SI-TH-99-3,hep-ph / yymmnn 
Skyrmions and Bags in the 2D-0(3) model 



G. Holzwarth* 

Fachbereich Physik, Universitdt Siegen, D-57068 Siegen, Germany 

(May 1999) 



Abstract 



Localized static solutions of the 2D-0(3) model are investigated in a rep- 
resentation with the 3-vector field split into the unit vector $ and the 
modulus <3?. As in the nonlinear version of the model this allows for the 
definition of a topological winding number B, and for the separation of 
the complete configuration space into distinct 5-sectors. For small values 
of the $ 4 -coupling strength the stable energy minima in these sectors are 
characterized by bag formation in the modulus field which in the standard 
cartesian representation of the linear 0(3) model would be unstable to- 
wards decay into the trivial B = vacuum. Stabilized by ^-conservation 
they exhibit a surprising variety of very appealing features for multiply 
charged systems. With the total charge bound into one common deep bag 
opposite ways of distributing the topological charge density inside the bag 
can be realized: Pointlike structures which retain the individuality of sin- 
gle constituents (or doubly charged pairs), or a deconfined charge density 
spread uniformly throughout the interior of the bag. It is suggested that 
this extension supplies a crucial link to overcome the unsatisfactory existing 
mismatch between multiskyrmion configurations and nuclear structure. 
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1. Introduction 



The linear 2D-0(3) model has been a favorite tool as a model field theory for a wide 
spectrum of physical systems. Commonly, the field vector is embedded into a eu- 
clidean manifold and parametrized in terms of three cartesian components (i=l,2,3). 
Topologically this manifold is trivially connected. For standard $ 4 potentials of the type 
— I) 2 the vacuum manifold is the 2-sphere S 2 . Addition of symmetry breakers can 
remove the degeneracy of the vacuum manifold. Imposing the constraint J2% &l = 1 de- 
fines the nonlinear 2.0-0(3) model. The homotopy S 2 — + S 2 then allows for a classification 
of static solutions by an integer topological winding index B which is the spatial integral 
over the time component of a conserved topological current. Nontrivial configurations 
with B ^ have found a direct application for the interpretation of the Quantum Hall 
Effect as charged excitations in two-dimensional spin systems 0. In the linear model, 
where configurations are not confined to the 2-sphere S 2 , the winding number B is not 
topologically protected. Therefore, beyond a critical magnitude for the coupling con- 
stants of symmetry breakers, static solutions with non vanishing winding 5^0 become 
unstable, and they unwind into the topologically trivial (B = 0) uniform vacuum. In 
time-dependent evolutions changes in B will frequently occur whenever the field passes 
through $ = at some point x, t in space and time. 

In certain applications of the model the length of the $- vector is an important degree of 
freedom such that the constraint to S 2 is too restrictive, but still, the winding number may 
correspond to a physically observable conserved charge. A prominent example is the chiral 
phase transition where the order parameter is the chiral meson field, its length measures 
the amount of spontaneous symmetry breaking, while the Skyrme-Witten conjecture || 
identifies the winding index with baryon number. Its conservation law should not be 
affected by the chiral transition. Then, in order to retain the model as a faithful image of 
the underlying physics, it is necessary to prevent unwinding of nontrivial configurations. 
To achieve this it is sufficient to exclude one poimt from the manifold on which the 
fields live. The natural choice is to remove the origin <& = 0, which is automatically 
accomplished in the radial-angular representation $ = because the angular fields 
are not defined at $ = 0. Then the manifold is nontrivially connected and the winding 
number is topologically protected. Thus, choosing the appropriate embedding for the field 
is an important part of the definition of the model. 

In lattice simulations the topological arguments based on continuity cannot be used 
because field configurations may differ arbitrarily between two neighbouring lattice points. 
However, the conservation law is easily implemented into the update-algorithm, by allow- 
ing only for configurations with specified value of B. Evidently, for that purpose it is 
necessary to be able to evaluate the winding number for each configuration, so it is nec- 
essary to work in the representation $ = $<!>, where the winding density is defined in 
terms of angles which may vary arbitrarily from one lattice point to its neighbours. 

So, in the following, we shall be dealing with all degrees of freedom of the linear 0(3) 
model, however in the form of a nonlinear 0(3) model supplemented by an additional 
modulus field $. In this framework it is possible to define and conserve the winding 
number of nontrivial configurations. Naturally, in a sector with specified B ^ 0, we will 
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find static solutions which do not exist in the topologically trivial cartesian embedding of 
the linear 0(3) model and which are characterized by the formation of a localized spatial 
bag in the modulus field. 

We shall discuss two versions of the model which differ in the definition of the cur- 
rent that enters the current-current coupling Skyrme term. The competition between a 
symmetry-breaking Zeeman term and the Skyrme term determines the spatial extent of 
the soliton solutions. Thus the two different versions will lead to a characteristic difference 
in the spatial structure of the resulting winding density distributions in the interior of the 
bags: A dominating Skyrme term deconfines the topological charge inside the bag while 
a dominating Zeeman term preserves the particle nature of individual charge units. 



2. The 2D-0(3) model with Skyrme- Zeeman stabilization 

We consider the 0(3) lagrangian density in 2 + 1 dimensions in terms of the dimen- 
sionless 3-component field <fr = with $ ■ $ = 1, 

C = F 2 Qa M #^* - A ($2 _ ff _ « (/o _ $ 3 ) _ ( a f) Pfl p^ . (1) 

Apart from the usual kinetic term this lagrangian contains the standard $ 4 potential for 
the modulus field $ to monitor the spontaneous symmetry breaking with (dimensionless) 
coupling parameter A, an explicitly symmetry-breaking ('Zeeman') coupling with (dimen- 
sionless) strength a, and a four-derivative ('Skyrme') current- current coupling p M p M for 
the conserved topological current 

(f = ■ {dA x d p *), (2) 

satisfying <9 M p M = 0. The parameter F 2 sets the overall energy scale, the length I may be 
absorbed into the space-time coordinates, so it sets the size of localized static solutions. 
In order to keep the uniform minimum of the potential for finite a at the vacuum value 
$ = fo we define 

f = fo ~ (3) 

then we have the a-independent boundary condition for the modulus $ at spatial infinity 
$ — ► /o for r — > oo. Of course, we are free to insert additional powers of the modulus field 
$ into the Skyrme and the Zeeman term, the above choice being motivated to minimize 
interference with the $ 4 spontaneous symmetry-breaking mechanism. This choice implies 
that as fo goes to zero (e.g. with increasing temperature) the typical size £ of static 
defects grows like f 1//4 . This may be physically not unreasonable (cf. e.g. the discussion 
in the 3-dimensional case in ||). (We shall discuss a different natural choice in the next 
section). Having fixed the ^-dependence of the lagrangian as given in (|IJ) and (0) we 
conveniently redefine the field and the parameters by 
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$ = $/ - 1 , F 2 = F 2 f , I=£fo 1/ \ (4) 



A = A/ 3/2 , a = a/- 3/2 , ~f = 1 - " (5) 

A 

omit the tildes in the following and absorb the ts into the length scale of space-time. 
Then we finally have 

£ = F 2 ^#0"* - ~ (V - 1 + - W - «(1 - $ 3 ) + ^ • (6) 

together with the a-independent boundary condition $ — > 1 for r — > 00. The constant 
ct 2 /(4A) has been added to set the ^-potential to zero at the minimum $ = 1. The 
nonlinear 0(3)-model emerges in the limit A —>■ 00 where $ is confined to the 2-sphere 
<& = 1. For this case the static solutions for the angular field 5> have been thoroughly 
investigated 0. They are characterized by fixed integer winding number B, and by a 
modulus field $(cc) which for very large A differs only minimally from its vacuum value 
$ = 1. We shall in the following denote them as S-skyrmions. 

For fixed a ^ and sufficiently small values of A the point $ = 1 is the only real 
minimum of the potential in @ , i.e. for A smaller than some critical value A c static 
solutions in the cartesian representation of the linear 0(3)-model collapse to the trivial 
vacuum $ = 1 . In a A- a diagram the line A c (a) separates two phase regions where 
for A > A c (a) static multi-skyrmions with winding number B can exist at local energy 
minima Eb, while for A < A c (a) only the global minimum <fr = 1 at E = survives. 
There will be different phase boundaries A c (a) for multi-skyrmions with different winding 
numbers B. 

In the <fr = representation of the linear 0(3) model the winding number is fixed 
(and can be held fixed in lattice simulations). This means that for A > A c (a) the local 
(multi-skyrmion) minimum Eb turns into a global minimum in the respective 5-sector 
and for A < A c (a) the collapse of the multi-skyrmions into the vacuum is prevented. 
Instead, the global minimum in a given 5-sector continues to exist for decreasing values 
of A. But we expect a qualitative change in the structure of stable static configurations 
of (Q) in the vicinity of the phase boundaries A c (a) leading to new types of solutions 
for A < A c which do not exist in the cartesian embedding of the linear 0(3)-model. 
The nature of this structural change is easily visualized: Destabilization and unwinding 
can only proceed through the field configuration passing through $ = at some point. 
Conservation of winding number therefore leads to formation of a spatial bag: In its 
interior the modulus field $ is close to zero, with rapid variation of the angular fields 
<E> such that the winding density is located inside the bag. Thus we may expect a close 
relation between the deviation of the scalar modulus field from its background value $ = 1 
and the density of the topological charge. We shall in the following denote these types of 
solutions as 5-bags. 
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3. Bag formation 



Within a given 5-sector the transition from skyrmions to bags is smooth. As A 
approaches A c from above creation of the bag begins through formation of a shallow 
depression in the modulus $ near the center of the skyrmions. As A passes the critical 
value the depression quickly deepens with the winding density accumulating inside. In 
fig. [I] this transition is illustrated in the B = 1 and B = 5 sectors. Plotted is the minimal 
energy per unit topological charge Eb(oc, \)/B as function of A for three different values 
of a. The full lines show the B = 1 sector. For a small value of a = 0.1 the transition 
near log A ~ 3.2 is most pronounced. Bag formation sets in rather sharply, with a sudden 
steep decrease in E B with decreasing A. For very large A above this critical value, E B rises 
very slowly towards a limit slightly above the Belavin-Polyakov (BP) monopole energy of 
47r [[J. For larger values of a the transition region gets smoothed out, still with a large 
gain in energy through formation of the bag. 




FIG. 1. Minimal energies per unit topological charge Eb{u,X) / B as function of log A for 
three different values of a = 0.1, 0.5, 1.0. Full and dashed lines show the cases B = 1 and B = 5, 
respectively. 

The dashed lines show the same features in the B = 5 sector. Comparing the energies 
for B = 1 and B = 5 for small A values near A ~ 1 which correspond to well-developped 
deep bags we see that E§ < 5Ei for all three values of a. This implies that the five 
topological charges are strongly bound in the compact B = 5 bag. It is only for the 
smallest value of a considered that with increasing A near log A ~ 1.75 the 5-bag breaks 
up into five individual 1-bags such that E 5 /5 closely follows E\ into the transition region 
where the bags disappear and the emerging five 1-skyrmions combine into pairs of two 
2-skyrmions plus one left-over 1-skyrmion (for log A > 3.2). For the larger a- values 
considered, E 5 /5 always stays below E\, so at no point there is a complete breakup into 
five individual 1-bags. 
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FIG. 2. Minimal energy Eb per unit topological charge B as function of B. (A = 10, a = 1). 

This holds for all values of B for sufficiently small A. Figure [| shows the energy per 
unit topological charge at A = 10 and a = 1. It is a monotonically decreasing function of 
B which converges towards approximately one half of the B = 1 energy; so for this A all 
B-bags are stable against decay into 1-bags. Altogether there emerges a phase diagram of 
remarkable richness. The skyrmion regions (A > 1) have been amply discussed, so here 
we only illustrate typical features of the region with well developped deep bags. (In order 
to get sufficient resolution the results are calculated for £ = 10 on an 80x80 mesh.) 




1 *m 




FIG. 3. 3D-view of the bags formed for B = 1 (a), and B = 20 (b), (for a = 1, A = 10). 



The 1-bag which is formed for a = 1 and A = 10 is shown in fig. ^a. Its total classical 
energy is E\ = 6.85. Already for this smallest value of B the bag has developped a 
basically flat bottom where the modulus field is very close to zero, with small numerical 
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FIG. 4. 3D-view of the angular field 0(x) = arccos $ 3 (cc) for S = 20 (A = 10, a = 1), 
fluctuating rapidly inside the bag of fig. [3|b around an average value of ir. 

fluctuations. To accommodate increasing numbers of topological charge B the bag size 
increases correspondingly, its depth being limited by $ > 0. As an example, fig.|3]b shows 
the case B = 20 for A = 10 and a — 1. The near-central cuts through the bag profiles 
plotted in fig. ^|a (for 5=1,2,4,6, ..,20) show that the surface thickness of the bags is 
basically independent of B, it is only the radius of the flat interior which adjusts to the 
increasing total charge. 

Throughout this deep bottom of the bag the angular fields display numerous strong 
local fluctuations. (Fig.f| shows the angle 9 = arccos $ 3 (for B = 20) fluctuating rapidly 
inside the bag around an average value of ir). Such fluctuations are energetically harmless 
because the large gradient terms are multiplied by the almost vanishing square of the 
modulus field. By this mechanism an essentially constant winding density po is achieved 
which extends over the whole flat interior of the bag (apart from small numerical fluc- 
tuations, see fig.^b). In the context of nuclear physics one would say that the nucleons 
have lost their identity inside the nucleus and have dissolved into a hadronic soup. In 
fig.^b the near-central cuts through the density profiles are shown which correspond to 
the same bags as given in fig.|]a. Again the surface thickness of the profiles is basically 
constant; it is, however, smaller than that of the bag profiles. So, for larger values of B a 
square slab closely approximates the density distributions. This is to be contrasted with 
the winding density of a B-skyrmion. Its angular fields closely resemble the BP monopole 
0(r) = 2 arctan(5/r) with winding density 

_ BO' sinG _ B 3 
PBP ~ A^r ~ vr(r 2 + 5 2 ) 2 ' ^ 

There is no central plateau, the center density decreases as p_ep(0) oc B~ x and the mean 
radius increases as R oc B. In contrast to this, fig.[5]b shows that for the B-bags the central 
density po quickly drops from its maximal value for B = 1 (p = 0.025 for A = 10, a = 1), 
then slowly converges with increasing B towards a constant value, which for this parameter 
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FIG. 5. Near-central cuts through the bag profiles (a), and the corresponding winding 
density profiles (b), (for B=l, 2,4,6, ..,20) for A = 10, a = 1. (Due to numerical fluctuations the 
positions of the bag centers may differ slightly for different values of B.) 



set is near po ~ 0.017. This -B-dependence reflects the interplay between the surface 
energy of the bag (which originates in the gradient terms of the modulus field) and volume 
energies which comprise the potential energy of the modulus field, Zeeman and Skyrme 
terms, and kinetic terms due to the gradients of the angular fields in the linear a-model 
part of the lagrangian (|) . 

To get an idea about the ('nuclear matter') density po m the limit B — > oo (where 
surface terms play no role) we could approximately replace the winding density defined 
in terms of the angular fields (|2|) by the deviation of the modulus field $ from its vacuum 
value 



P « Po(l - <*>) (8) 

and obtain the bag radius R from TrR 2 p Q = B. The lagrangian (^ then shows that 
creation of a deep bag ($ ~ inside, $=1 outside) requires a volume potential energy 
density of 

A + 2a „ 
ev = + 1 ap v (9) 

Ignoring kinetic contributions for the moment, for fixed B the interplay between these 
two volume terms determines the average density po as 




. \ + 2a , , 

which for A = 10 and a = 1 is p = V^/f = 0.0173. 

The additional pressure due to the kinetic volume terms will further lower the central 
density, but comparison with the numerical result in fig.^jb shows that their influence 



8 



must be small. This is what we might expect due to the smallness of $ inside the bag, 
although very large gradients in $ could compensate for it and cause a noticeable lowering 
of po, especially for small values of a. Maximal importance of kinetic terms would occur 
if the hadronic soup inside the bag would consist of fermions. Then, in Thomas-Fermi 
approximation for the density p — Pf/(4tt) we would have for the energy per particle 

6 ^ ^~ | ev (11) 



Po 3 po 
Variation with respect to po then leads to 

^ „ 3/2 A + 2a 2 2 



'0 



U 2 



+ tapl = 0. (12) 



For A = 10 and a = 1 the central density thereby is lowered from 0.0173 to p = 0.0120. 
Comparison with fig||3 apparently rules out the conjecture of a fermionic soup. (For small 
values of a the effect of the Thomas-Fermi term becomes stronger: for A = 10, a = 0.1, 
£ = 10, it lowers the result of ( |l0l) po = 0.050 to po = 0.016, while the numerically 
obtained density for B = 20 is near 0.04). 



4. Individual particles inside the bag 

An interesting alternative way to define the Skyrme term in the lagrangian is to replace 
in ([I]) the topological current p M (|2|) by the corresponding form in terms of the full vector 
field $ 

T» = — ■ {d v <S> x d p &) = $V • (13) 

This introduces a factor $ 6 into the Skyrme term. Due to this high power of $ there 
again is very little interference with the $ 4 spontaneous symmetry-breaking mechanism 
of the ^-potential. However, apart from the fact that the length £ now scales with / 
as fo tends to zero, this alternative definition leads to a situation where, for fixed /o, in 
spatial regions with small $ (i.e. in the interior of bags) the Zeeman term dominates 
the Skyrme term. This results in the stabilization of pointlike 'particles' inside the bag. 
Speaking again in the language of nuclear physics, one would say that the nucleons retain 
their individuality inside a bag which binds them together into a common 5-nucleus. 

Choosing as before the parameters A = 10, a = 1 which favor the formation of deep 
bags, we find the energy minimum for B = 1 at E\ = 1.152. In the B = 2 sector 
we find a minimum near E 2 /2 = 0.943 in which both charges share one common bag, 
but the individual charges stay apart from each other in two separate pockets. This 
'l + l'-configuraton is shown on the left side of fig|| There is, however, in this sector a 
lower minimum at E 2 /2 = 0.663 where the bag consists of only one single pocket which 
houses one doubly charged structure, i.e. the model favors the formation of pairs. This 
'2'-configuration is shown in the center plot of fig.|6]. For B = 3 we find the expected 
'l+2'-configuration at E 3 /3 = 0.697, a pair and one unit charge tightly bound in two 
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FIG. 6. 3D-view of the bags and corresponding winding densities for the '1+1' and '2' 
configurations observed in the B = 2 sector (left and center) and the '1+2' configuration in the 
B = 3 sector (right), if the Skyrme term in (Q) is multiplied by $ 6 as follows from the alternative 
definition (13). (a = 1, A = 10). The energies per unit charge are Eb/B = 0.943,0.663,0.697, 
respectively. 



separate pockets of a common 3-bag, as shown on the right of fig.|| We do not find a 
minimal-energy ^'-configuration with a threefold charge sitting in one single pocket. The 
same situation is observed in the B = 4 sector where again we do not find the single 
pocket with fourfold charge, but instead the '2+2' double pair separated in two pockets 
bound in one common 4-bag at E^/A = 0.543. So, it appears that the situation resembles 
very closely the pair formation as it was observed M for skyrmions for A — > oo. As in that 
case, the '1'- and ^'-configurations serve as building blocks which here are bound together 
in one deep bag. Due to this energetically favored pair formation the resulting 'nuclei' will 
be characterized by odd-even staggering of their binding energies for higher values of B. 
Practically, of course, the actual configuration which is finally reached through numerical 
relaxation depends on the starting set and the sequence of pseudo temperatures in the 
Metropolis cooling algorithm. As an example we present in fig.|7| a result for B = 20 
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FIG. 7. 3D-view of the bag formed for B = 20, and the corresponding winding density p, 
if the Skyrme term in ([EJ) is multiplied by <J? 6 as it follows from the alternative definition fll3|). 
The parameters (a = 1 and A = 10) are the same as in fig. |3|b. 

which shows a highly deformed bag with 15 pockets which houses a 10 x' 1' + 5 x' 2' 
configuration. The calculation is done for the same parameter set as used in fig.^b, so 
a direct comparison shows nicely the contrast between pointlike particles and deconfined 
charge in the interior of one big bag. 

5. Conclusion 

We have investigated here the 2.0-0(3) model in a representation where the 3- vector 
field $ is split into the unit vector <3E> and the modulus <3>. This allows for the definition 
of a topological winding number B, and for the separation of the complete configuration 
space into distinct 5-sectors. For small values of the $ 4 -coupling strength A the sta- 
ble energy minima in these sectors are characterized by bag formation in the modulus 
field. In the standard cartesian representation of the linear 0(3) model such configu- 
rations would be unstable towards decay into the trivial B = vacuum. Stabilized by 
-B-conservation they exhibit a surprising variety of very appealing physics for multiply 
charged systems. For decreasing A multi-skyrmions get bound into one common deep bag 
like nucleons get bound into one nucleus. Depending on the competition between Skyrme 
and Zeeman energy two opposite ways of distributing the topological charge inside the 
bag can be realized: Pointlike structures which keep the individuality of single nucleons 
(or doubly charged pairs) inside the nucleus, or a deconfined charge density spread uni- 
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formly throughout the interior of the bag. This latter case suggests a very close relation 
between the charge density p and the modulus field $ (which for sufficiently large B is 
p = po(l — $)) and, correspondingly, an effective description through a density functional 
for p. This is a remarkable possibility because it gets rid of the angular fields altogether 
which form the basis for the definition of p ! It may be understood by the fact that in 
these configurations the angular fields are very rapidly fluctuating functions and a spatial 
coarse graining procedure will eliminate them. This is reminiscent of the fact that mean 
field or density functional methods for nucleons in nuclear physics work very well without 
keeping explicitly the dynamical pionic degrees of freedom. Similarly, in a recent analy- 
sis of soliton formation in the Nambu-Jona Lasinio model || it was found that without 
enforcing the chiral circle condition stable minima are characterized by vanishing pion 
field. 

Of course, it will be most interesting to extend the present investigations to the 3D- 
0(4) model, where a lot of effort has gone into exploring 3D-SU(Nf) multi skyrmions |7J 
and their possible relevance for the structure of nuclei. Their spatial structure (which 
for large B looks somehow like buckey balls) has little in common with our naive picture 
of a nucleus. On the other hand it has been known for a long time || that inclusion 
of a scalar a field is important for the attractive part of the skyrmion-skyrmion force. 
We should also rather expect a mechanism as described in sect. 4 where the nucleons 
keep their individuality inside the bag. For dimensional reasons in the 3.0-0(4) model 
there is no need for the Zeeman-type coupling to monitor the scale of the structures; it 
emerges directly from the competition between the second order (non linear a) term and 
the (fourth order) Skyrme term. Dominance of the Skyrme term (which may arise in 
the interior of the bag due to the factor <3> 2 multiplying the nonlinear cr-term, or due to 
increasing temperature which lowers the coefficient of the whole second order term) will 
again lead to deconfined baryon density inside the nucleus (or nuclear matter in the limit 
B — > oo). In any case, our present results strongly suggest that for a description of nuclei 
in terms of multiply charged solitons in chiral meson fields the constraint $ = 1 is too 
restrictive. 

Creation of spatial regions with disoriented chiral condensate (DCC) in the course of 
a chiral symmetry-breaking transition has been studied in the framework of the linear 
0(4) model in trivial topology ||. It will be of interest to investigate how the existence 
of the various phase regions explored here for the simple 2.0-0(3) model may affect 
the conclusions to be drawn for defect and DCC formation in case of the chiral phase 
transition [|K] . 



Finally we should note that all results reported here have been obtained through 
numerical relaxation of a lattice functional through some Metropolis cooling algorithm. 
This naturally poses the question about the continuum limit of these results. Apparently, 
(as is directly obvious from fig.|4]), this is not a trivial matter because B conservation 
relies on the removal of one single point ($ = 0) from the manifold on which the model is 
defined. But this is the usual mathematical difficulty encountered in the transition from 
a granular to a continuous density distribution. And, again as usual, we have no proof 
that we really have obtained the lowest energy minima in the respective O-sectors. 
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